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Abstract 

For a class of systems of nonlinear and nonlocal balance laws in several space dimensions, we 
prove the local in time existence of solutions and their continuous dependence on the initial 
datum. The choice of this class is motivated by a new model devoted to the description of 
a metal plate being cut by a laser beam. Using realistic parameters, solutions to this model 
obtained through numerical integrations meet qualitative properties of real cuts. Moreover, 
the class of equations considered comprises a model describing the dynamics of solid particles 
along a conveyor belt. 
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1 Introduction 

We are concerned with a system of n balance laws in several space dimensions of the type 

J dtUi + AYVxVi{t,x,Ui,-d *u) = *u) ._ . , 

I Ui{Q,x) =Ui[x) i- 

Here, t G [0, +oo[ is time, x € is the space coordinate and u = (mi, ..., tt„), with Ui = Ui(t, x), 

is the unknown. The function r? is a smooth function defined in R'^ attaining as values m x n 
matrices, so that 

C2(R^;R’”^”), {-&*uit)){x)= [ - p u(<, ^ , (r? * u(t)) (x) G R™ . 

The flow ip = with ipi{t, x,Ui, A) G R^, and the source $ = with 

X, Ui, A) G R, have the peculiar property that the equations are coupled only through the 
nonlocal convolution term * u. 

The driving example for our considering the class (ED is a new model for the cutting of metal 
plates by means of a laser beam, presented in Section[3] A sort of pattern formation phenomenon, 
typical of various nonlocal equations [7], accounts for the formation of the well known ripples 
whose insurgence deeply affects the quality of the cuts. In fact, two type of lasers are mainly 
used in the cutting of metals: CO 2 lasers and fiber lasers. The former ones are more powerful and 
more precise, but also more expensive. Recent technological improvements are apparently going to 
allow also to the cheaper devices of the latter type to cut thick plates, nowadays treated typically 
with CO 2 lasers. Unfortunately, a typical drawback of fiber lasers is that along the cut ripples 
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are generated, see [laiiEiiss]. The modeling of these ripples often relies on the introduction of 
imperfections in the metal or of inaccuracies in the laser management, see also [251122]. Here, 
using realistic numeric parameters, we obtain the formation of a geometry similar to the ripples 
observed in industrial cuts. We remark that in the present construction neither the initial data 
nor the parameters in the equations contain any oscillating term. 

Furthermore, in Section 01 we slightly extend the model introduced in m to describe the 
dynamics of bolts along a conveyor belt. The resulting equations fit in the present framework and 
is proved to be well posed. 

Besides, we also note that several crowd dynamics models considered in the literature ht 
into (EH), e.g. EiniiiiKiH]. 

The particular structure of EH) allows to prove its well posedness. Indeed, for small times, 
system EH admits a unique solution u = u{t,x). Moreover, u is proved to be a continuous 
function of time with respect to the topology and an L^-Lipschitz continuous function of the 
initial datum u. In all this, the particular coupling among the equations in EH plays a key role. 
At present, the well posedness of general systems of balance laws in several space dimensions is a 
formidable open problem. In the present work, the functional setting is provided by L^nL°°nBV, 
as usual in the framework of nonlocal conservation laws. The existence result is obtained through 
a careful use of the general estimates [TOl [21] . They provide the necessary analytic tool to apply 
Banach Contraction Theorem. 

A preliminary result related to Theorem 12.21 below is presented for instance in [T] , see also [3] . 
There, the existence of solution to EH in the case $ = 0 is obtained proving the convergence (up 
to a subsequence) of a Lax-Friedrichs type approximate solutions. Note however that differently 
from the present situation, in the case considered in [1] , positive initial data yield positive solutions 
so that the norm is conserved. 

We remark that most of the results related to nonlocal balance laws are currently devoted to 
conservation laws, i.e., to equations that lack any source term. Here, we allow for the presence 
of source terms that can be nonlinear in both the unknown variable u and the convolution term 
'd * u. The unavoidable cost of this extension is a local in time existence result, as shown by an 
example in Section |2l 

Nonlocal conservation and balance laws are currently widely considered in various modeling 
frameworks. Besides those of crowd dynamics, laser cutting and conveyor belt dynamics considered 
above, we recall for instance granular materials, see |5], and vehicular traffic, see [1]. For a different 
approach, based on measure valued balance laws, we refer to m- 

The paper is organized as follows: the next section is devoted to the analytic results. Section[3] 
presents the laser cutting model, its well posedness and some qualitative properties with the help 
of numerical integrations. Conveyor belts dynamics is the subject of Section01 All analytic proofs 
are postponed to the last Section 0 

2 Analytic Results 

Throughout, we denote by grad^. /, respectively diVa, /, the gradient, respectively the divergence, 
of / with respect to the x variable, with x G . All norms in function spaces are denoted with a 
subscript indicating the space, as for instance in ||w(t)||Li(R]v.R,.)- When no space is indicated, the 

norm is the usual Euclidean norm in for a suitable k, as for instance in \\u{t, a:) ||. Throughout, 
we fix the non trivial time interval I = [0, T], For any U > 0, we also denote Ujj = [—U, U]. 

Our starting point is the definition of solution to EH, which extends [H Definition 2.1] to the 
case of balance laws. 

Definition 2.1. Fix a positive T. Let u G A map u : [0,T] —L°°(]R^,]R") is a 

solution on [0, T] to (11.11) with initial datum u if, for i = 1,... ,n, setting for all w gR 

(pi{t, X, w) = ipi (t, X, w, (i9 * u)(t, x)) and ^i(t, x, w) = (t, x, w, (d * u){t, x)) 
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the map u is a Kruzkov solution to the system 


dtUi + divj; ipi{t, X, Ui) = X, Ui) 

Ui{Q,x) = Ui{x) 


( 2 . 1 ) 


Above, for the definition of Kruzkov solution we refer to the original EOl Definition 1]. 

We are now ready to state the main result of the present paper. 

Theorem 2.2. Assume that there exists a function A S (C° nL^)(/ x x ]R+;R+) such that: 

(if) For anyU>0,(pG (C^ n x x Uu x and for all t G T, x G , 

uGUu, A GUff 



\\gi&d,^ip{t,x,u,A)\\ , 

' 

diva; (p(t, a;, M, A) , 

max < 

grad,, A\Y,,ip{t,x,u,A)\\ , 

grad„ grad^ 'pft, a;, u. A) , 

1 

grad^ V5(t, a;, u, A) , 

grad^ ip{t, X, u, A) 


($) For any U > 0, ^ G (C^ n x x Wy x R") and for alltGl,xG R^, uGUu, 

AGUff} 

max I ||<i)(t, a;, u, A)|| , Hgrad^, $(t, x, u, A)|| | < X{t,x, U). 

(i9) i? e C2(R^;R™X’"). 

Then, for any positive C there exists a positive T^, G I and positive C, C such that for any datum 

UG (L^nL°°nBV)(R^;R”) with ||tti|lLi(R«;B") ^ ll^i|lL~(R«;R") TV(Ui) < C, (2.2) 

problem CH) admits a unique solution 

uGC° ([0,T*];L^(R^;R”)) 

in the sense of Definition \2.1[ satisfying the bounds 

lhW|lLi(R«;R") lhW|lL“(R^;R") TV(u(t)) <C, 

for all t G [0,T*]. Moreover, if also w satisfies (12.211 and w is the corresponding solution to (HI, 
the following Lipschitz estimate holds: 

||u(t) — w{t) if kl, I I'M — '*fi|lLl(R«' iR") • 


The proof is deferred to Section O Observe that the whole construction in the present paper can 
be easily extended substituting the convolution id * u with a nonlocal operator having suitable 
properties that comprise those of the convolution, as was done for instance in mm- 


A natural question arises, namely whether the above result can be extended to ensure the 
global in time existence of solutions. In this connection, consider the following particular case 
of (dH) 


dtu = (u* r])u 
u(0, a;) = 1. 


Here, n = 1 and m = 1 while N does not play any particular role. Moreover, ij G C^(R^;R) is 
non negative and satisfies ? 7 (a:) da; = 1. The solution is u(t,x) = 1/(1 — t), which exists only 
up to time t = 1. The above example (lOl) admits an explicit solution but does not fit into the 
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setting of Theorem 12.21 since the initial datum is not in On the other hand, setting 

iV = 1, the similar problem 


dfU = {u* r])u 
u{0,x) = i^{x) 


where 


tp{x) 


1 |x| G [0,1] 

1)3)4 1^1 ^ 

0 |a:| G [2,+oo[ 



apparently has a qualitatively analogous blow up pattern, as shown by the numerical integration 
displayed in Figure |2j To obtain it, we use an explicit forward Euler method, with space mesh 
Aa; = 10“^ and time mesh At = 10“^ on the space domain [—3, 3] and for t G [0, 1.05]. The 



Figure 1: Numerical integration of (12.41) for t £ [0, 1.05]. The values of the 
norm of the solution is plotted vs. time in Figure [2] 

graph of the norm of the numerical solution to (12.41) is in Figured It is straightforward to see 
that dZl fits into the framework of Theorem 12.21 setting 

N = \ 

^ ^ I ‘p{t,x,u,A) = 0 

, ^(t.x.u.A) = ih(x)uA. 

m = 1 

The requirements {(p) and ($) are easily seen to be satisfied. 

3 A Laser Beam Cutting a Metal Plate 

A thin horizontal metal plate can be cut by means of a moving vertical laser beam. More precisely, 
the laser energy melts the metal along a prescribed trajectory. A wind, suitably provoked around 
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Integral of u 



t 


Figure 2: norm of the solution to (12.41) . suggesting a blow up at finite time, 

similar to the solution to (E31) . 


the beam, pushes the melted material downwards. For its industrial interest, this phenomenon 
is widely considered in the specialized literature, see [T31 [TH [13 [13 [23 HZl [23 IS21 133], while 
information specific to the cut of aluminum are for instance in m- A phenomenological description 
of the whole process can be summarized as follows. We fix a 3D geometric framework, with the 
laser beam parallel to the vertical z axis, see Figure[3 left. The trajectory of the laser is prescribed 
by the map xl = XL{t). We distinguish the height hs of the solid metal and that of the melted 




Figure 3: Left, reference frame with respect to the metal plate being cut. The 
laser beam is parallel to the z axis, while the plate lies on the z = 0 plain. 

Right, the distinction between the melted part hm. and the solid one ha. 

part, denoted hm, see Figure [3 right. 

A ID system of balance laws is used to describe the dynamics of the melted and of the solid 
material in |5| . Flere, we present a description of this dynamics by means a 2D system of balance 
laws of the form: 

I ^thm T div3;(hm V^) = £. Ij 

1 dths = -C. ^ ’ 

The vector V = V(t, x) describes the projection of the melted material velocity on the horizontal 
(a:, ?/)-plane. Its modulus must depend on the wind speed w = w(t,x), which is centered at the 
laser beam sited at x = XL{t). Its direction depends on the geometry of the melted metal and of 
the solid surface z = H{t,x), where H = hg + hm- The source term C is directly related to the 
laser position and intensity: it describes the net rate at which the solid part turns into melted. 
Also C depends on the metal geometry, since the heat absorption is strictly related to the incidence 
angle between the moving melted metal surface and the vertical laser beam, see Figure [dj left. 
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Here, we posit the following assumptions: 


V = {w{t,x) - Tghjn) - 

1+ ||grad^(7?*iJ) 

The term with the coefficient Tg in (E2I) is related to the shear stress, inspired by [SI [33] ■ The 
denominator in (13.21) is due to a (smooth) normalization of the direction — grad,j,(r? * H) oi the 
average steepest descent along the surface z = Indeed, the convolution kernel t] is chosen 

smooth, compactly supported and with total mass 1, so that grad^, [r] * H{t)) (x) is the average 
gradient at position x and time t of the surface z = Hlt,x). 

In (13.31) . the numerator i = i(t, x) is related to the laser intensity. It can be reasonably described 
through a compactly supported bell shaped function centered at the location of the moving focus 
of the laser beam. The denominator is the squared cosine of an averaged incidence angle of the 


-grad,,(??*g) 
Jl + ||grad^(7? + iJ)| 


(3.2) 

(3.3) 




Figure 4: Left, the incidence angle a of the laser beam on the surface z = 
H(t,x). Right, a possible profile for the functions W and T in (13.41) . 


laser on the surface z = H{t,x), see FigurejU left. In fact, 


/ 


COS a = 


[-d^,H -d^,H 1 ] 


\[-d,,H -d,,H 1] [0 0 1] 


1 + l|grad,^^ir 


V 

For the wind function w = w{t,x) and for the laser intensity function i = i(t,x) we choose a 
dependence on the form 


w{t,x) =W (^\\x — XLit)\\j and i{t, x) = I (^\\x — XLit)\\j (3.4) 

where both maps W and I have the form in Figure [d] right. More precisely, in the real setting 
under consideration, the diameter of the support of W is a few times larger than that of I. 

We stress that the present model describes how the laser beam digs a block of metal along its 
movement, i.e., it describes the dynamics of the melted metal and the profile of the solid material 
during the passing of the laser beam. At the physical level, the actual formation of the hole makes 
the melted material fall and, essentially, disappear. At the analytic level, the appearance of the 
hole causes major discontinuities that can hardly be described within a model of the form (13.11) . 
Therefore, we provide with an initial datum 


hs{0,x)=h° and hm{x) = 0 


(3.5) 
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where the constant /i^ is the uniform thickness of the plate under consideration. Then, we interpret 
the region where hs{t,x) < 0 as the region where the cut is accomplished. 

As a result we obtain the following model: 


- grud^jr] * H) 

1 + ||grad,,(77*i7) 


1+ ||grad,^(??*7J) 

U = hg T 

\ 

To apply Theorem 12.21 to the model (13.61) . a formal modification is necessary. Indeed, we 
introduce a cutoff function 

rgit,x) =TgS (\\x-XLit)\\^ where 5(e) = | J ^ ^ +oo[ 

for a smooth S and suitable (large) r and i?, with r < R. We thus obtain 


i{t, x) 


1 + ||grad ,^(?7 * H)\\ 

(3.6) 


dthm + diva; 


^w(t^x)h^ Tg(h^) 
i(t, x) 


dthm + diVa 
dths = -- 


- {w(t,x)hm - Tg{t,x)(hm)‘^) grada,(r 7 * H) 


i(t, x) 


V^l+ ||grad,,(7?*iJ)| 


i(t, x) 


1+ ||grada,(77*^^)| 


(3.8) 


1 + ||grada,(?7 * H)\ 

H — hg T hjYi 


When used with real data, the two problems (13.6p and (13.81) are indistinguishable. 
Proposition 3.1. The model (13. 8|) fits into (HID setting: 


N = 2 
n = 2 
m = 2 

U\ = hjYi 

U2 = hg 


dx^ri(x) d^^r](x) 
dx2V(x) dx^vix) 


ipi(t,x,ui,A) =(w(t,x) -Tg(t,x)ui) 

ip 2 {t,x,U 2 ,A) = 0 
$i(t,a;,u,A) = --^=A=i(t,x) 

^ 2 it,x,u,A) = --^=A=i(t,x), 


where w,i are defined in (1^ and Tg in dSlD- Moreover, if 


XL e (C^nW2’“)([0,f];]R2)^ >V,I,5 gC^ 


and rj G cl 


(3.9) 


for a positive T, then, assumptions (<p), (^), ("&) hold. 


The proof is deferred to Section [5] The above Proposition 13.11 allows to apply Theorem 12.21 to 
model dUD, ensuring its well posedness. 


3.1 Numerical Integration 

The model (13.61) . fed with realistic values of the various parameters, is able to reproduce the rising 
of ripples. The following numerical integrations show this qualitative feature. 

We use below the numerical method presented in [1], where it is proved to be convergent up 
to a subsequence in the case of a system of nonlocal conservation laws. As it is usual, we deal 
with the source terms by means of the fractional step method, see for instance [HI Section 12.1]. 
In other words, we use a Lax-Friedrichs type algorithm for the convective part and a first order 
explicit forward Euler method for the ordinary differential equations arising from the source terms. 
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The computational domain is the rectangle [0, 40] x [—2, 2], entirely contained in the metal 
plate to be cut (all lengths being measured in millimeters). The mesh size is 5 • 10“^ along both 
axis. The integration is computed for t G [0, 1], time being measured in seconds. 

The laser trajectory is 


XL{t) 


(3, 0) t G [0, 0.1] 

(3 + 40(t-0.1), 0) t e]0.1, 1] 



meaning that for t G [0, 0.1] the initial hole is drilled centered at (3, 0), in the interior of the metal 


Figure 5: Numerical integration of (13.61) with the data and parameters provided in S I3.ll see [511^. 
These are the contour plots of the solid metal level hs in the interval [0, 4.5] over the domain 
[0, 40] X [—2, 2] millimeters. The inner circle (appearing as an ellipse due to the different scales 
on the two axis) is the support of the laser beam. The outer one is the support of the wind. At 
time t = 0.1, the initial hole is terminated and the laser beam starts moving rightwards. Note the 
formation of “ripples”, i.e., the sides of the cut are not flat but present an apparently regularly 
oscillating profile. Neither data nor parameters are “pulsating”. 


The wind and laser functions are given by (15a setting 


W(0 


2(l-(^)) for lien <1.2 


corresponding to a laser beam with radius 1.2 mm, see m- The radius of the surface where the 
wind blows downward is 3 times that of the laser beam. Besides, we set Tg = 4. The convolution 
kernel is 


r]{x) 


1 

/r2 viy) dy 


fl{x) 


where fi{x) 



As initial datum, we choose 


for ||x|| < 2.4 . 


hm{x) = 0 and 


hs{x) = 4.5 for all x S , 


representing a flat metal plate 4.5 mm thick, see m- 
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t=1.00 



Figure 6: Left, a real piece of metal with a hole and a cut made by a laser 
beam. Right, the result of the numerical integration of the model (I331l as in 
Figure [S] but plotted with the same scales along the two axis. 

The result of this integration is in Figure [SI which displays the contour plot of the surface z = 
hs{t,x) remaining after the cut, restricted to the interval [0, 4.5]mm. The white part corresponds 
to a level below 0 and should be understood as expelled, corresponding to the cut. Remark the 
oscillations arisen along the sides of the cut. No parameter and no datum in the integration 
oscillates, nevertheless, the solution displays these sort of “ripples”. 

4 Materials Flowing on a Conveyor Belt 

A macroscopic model for the flow of materials along a conveyor belt is presented in [12]. The 
material consists of a large number of solid identical particles, called cargo. From a macroscopic 
point of view, the cargo state is identified by a density p = p{t, x), where t is time and x = (xi, X 2 ) 
is the coordinate along the conveyor belt. The industrial interest behind these modeling efforts 
is motivated by the need of an efficient management of specific parts of the production process. 
A standard example is the pouring of newly produced bolts in boxes. In this case, a selector is 
positioned on the belt to drive the bolts in a short segment of the belt, so that at the end of the 

conveyor they fall in their boxes, see [T^j and figureS] left. For other references on these modeling 

issues, both from the microscopic and macroscopic points of view, see for instance [l^j related to 
conveyor belts in mines, or [niiis] and the review m- 

With the notation in [H Section 3], a macroscopic description for the cargo dynamics is 
provided by the equation 

dtP + dYVa:{^p{y^*'°‘^{x) ^ H{p- Prnax)'I-{fi)^ =0. (4.1) 

Here, is the time independent velocity of the underlying conveyor belt. The fixed positive 

Pmax is the maximal cargo density and H is the usual Heaviside function. The term 

= (4.2) 

^1+ ||grad,,(? 7 *p)f 

describes how the cargo velocity is modified when the maximal density is reached: particles move 
towards regions with lower average cargo density, rj being a positive function with integral 1 , 
so that 77 * p is an average cargo density. Further details are available in [m Section 3], where (j4.1|) 
is supplied with suitable boundary conditions along the sides of the conveyor belt. The numerical 
study therein shows a good agreement between the solutions to (ED and real data. 

Next, we slightly modify 6ID- The conveyor belt is described by the strip |a: 2 | < i. First, we 
replace the Heaviside function by a regularization 

H>^ e C2(IR; [0,1]) with = H{^) with |^| > p. (4.3) 
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Then, we modify so that it incorporates the upper and lower conveyor boundaries. To 

this aim, we introduce the vector field b{x) S C^(IR^;R^), see Figure El right, such that: 



Figure 7: Left, a conveyor belt with a selector restricting the possible path of 
the carried cargo. Right, geometry and notation of the conveyor belt. 


{b{xi,£))^ =—e Va;ie]R, (b(a;))^ = 0 Va; G , 

(b{xi,—£))^= e Va;iGR, (b(x))^ = 0 Va; G with |a ;2 — ^| > ^ or |a ;2 + ^| > <5 . 

We therefore obtain the equation 

dtp + diva; ^p (^u®*“*(a:) + b(x) + H^^{p - Pmax)I{p)^'^ = 0 (4.4) 


which describes the cargo dynamics along the conveyor belt. Thanks to the framework provided 
by Theorem 12.21 we can incorporate in the model also the cargo source and sink. Indeed, we 
assume that the solid particles are poured on the belt in a region, say, = [0, a] x [—£, £\ and fall 
out of the belt in the region Rout = [T — a, L] x [—£, £]■ To this aim, for a positive T, we introduce 
the source and sink functions 


4'„GC2([0,f]xi?„;R+) 


with spt •) C Rin for all t G [0, T] 


4'oiit G C^dO,?] X X ]R;R+) with 


spt '^out{t,-,p) C Rout for all t G [0, T] and p G M 
4'o«t(-,-,p) = 0 for all p < 0 


(4.5) 


The function is the rate at which particles are poured in Rin, while '^out describes the outflow 
from the belt. 

We can assume that the belt is initially empty, thus we obtain the following Cauchy Problem, 
where we set v = + b, 


{ dtp + dlYa: (^p {v{x) + H>^{p - Pmax) I{p))'^ = x) - out(t, X, p) g. 

p{t,0) = 0 

Proposition 4.1. Fix positive T,£, L, p,, pmaxT£,£ with e > e. Let B = [0,L] x [—£,£] be the 
conveyor belt. If v G C^(R.^;K^) is such that 

{v{0,x2))^ > 0 yx2 €[-£,£] {v{L,x2))^ < 0 yx2 €[-£,£] , . 

{v{xi,-£)).^ > £ Va:iG[0,L] (v{xi,£)).^ > £ Vxi G [0, L] , 

is as in (ig , X is as in (14.21) and d/in, d/out are as in (ig . then there exists a positive T* 
such that problem (EH) admits a solution on the time interval [0,T*]. Moreover, this solution is 
supported in B for all t G [0, T*]. 

The proof is deferred at the end Section El 
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5 Technical Details 


The proof of Theorem 12.21 consists of several steps. We briefly describe here the overall formal 
structure. 

Fix a positive T G / and let / = [0,T]. Introduce the map 

T : {w, u) ^ u 

where u = (ui,..., Un) and its Tth component ui solves the nonlinear balance law 


dtUi + divj; * w) = ^i{t, x,Ui,'d * w) 

Ui{0,x) = Ui{x) 


(5.1) 


for i = 1,... ,n. By construction, solving dm is equivalent to solving the fixed point problem 
u = T{u,u). The core of the proof thus consists in choosing T and suitable subsets 

WC C° and UcL^(R^;R"), 

see (ED), so that 


(i) 

y{w, ft) G w X u, 

T(w,u) is well defined. 

(ii) 

y{w, fi) G w X u, 

T{w, ft) is in W, 

(iii) 

VftG u, 

w —>■ T(w, ft) is a contraction. 

(iv) 

Vw G W, 

ft —>■ T (w, ft) is Lipschitz continuous. 

(v) 

V(w;, ft) G W X U, 

t —>■ (T('ut, ft)) (t) is continuous. 


Steps 1 and 2 in the proof below give (i). The a priori bounds proved in steps 3, 4, 5 and 6 
ensure (ii). The key estimate (I5.19L which has the form 

||T(w , u) — T(ui )W)||po(/.Li(RiV.Rn)) < C>(1) ^11'^ ~ IIc0(/;L1(R«‘;R»)) 

and is proved in Step 7, shows that (iii) holds for T small. The statement (iv) is obtained in 
Step 9 through an estimate of the form 

\\T{w,U ) — T{w,u )||co( 7 ;Li(R^;R")) — 11^ IIli(R^;R")> 

see (15.201) . Finally, (v) is the content of Step 4, see (I5.9I) - (I5.10L used also in the proof of (ii). 

Once the statements (i), ..., (v) are obtained, the proof of Theorem l2.2l is essentially completed. 
Proof of Theorem 12.21 Throughout, we use the standard properties of the convolution product 
and, in particular, the following bounds. If d satisfies {'&) and u G L^(R'^;R"), then 

ll'i? * lt||Loo(7xR«‘;R") — II'*^IIl“(R^;R’"X") II'*^IIc'’(/;L1(R^;R)) 

which is a straightforward generalization, for instance, of [5J Theorem IV. 15]. By (i?), without 
any loss of generality, we may assume that 

II^v|Ili(R~;R) - for all j = and i = l,...,n. 

This requirement simplifies several estimates below, since it ensures that 

Ui{x) G Uu for alH = 1,..., n and x G R^ => (z? * u){x) G for all x G R^ . 
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1: Notation and Definition of T. Fix positive iF, U, U, R and R with 

II^IIl’-(r''^;R'‘) ^ -R < -R I II'^IIl°°(r'^;R’*) ^ U < U and TV(u) < K. 

Introduce the closed sphere centered at the initial datum u with radius R and its intersection 
with BV as follows: 

Bi^i{u,R,U) = e L^(]R^;R"): ||m - ullLi(Riv.Rn) <-R and u(a;) G Wyj 
RLinBv(u,R,R,R) = {uGB-^iiu,R,U): TYiu)<K}. 

For any positive T G I, denote I = [0, T] and define the map 

T : CO(/;RLi(h,R,t/)) x R,U, K) ^ C° {l; B^iiu, R,U)) 

w , u u \ / 

where the function u = {ui, ..., is such that for i = 1,..., n, Ui solves (|5.1I) . We equip the 
Banach space C° (/; L^(K^;K”)) with its natural norm 

ll^llco(/;Li(R'^;H")) = II^(^)|Il1(R'^;R") ’ 

and the metric space Bi,ir\u\'{u,r, K) with the L^~distance. Denote below 
VLt = / X X R and = / x x Uu . 


Moreover, we set 

A{t,U) = ||^(•) •) R)||Ll([ 0 _t]xRW;R) (5-3) 

so that A(-, U) G C°(J; K) is non decreasing, bounded and A(0, U) = 0 for all U G K'*'. 

Throughout, we denote by C a quantity dependent only on A and on the norms in (c^), ($) 
and (d), but independent of T, R, U, R, U and K. Similarly, Cjj is a constant depending only 

ll‘/5|lw2.“(/xR"xWt/xW5‘;R"X™) II *^11(/xR''xM[/xWJf ;R")’ 

2: Problem (15.11) Admits a Solution. Note that if i 7 ^ j, equation (15.111 is decoupled from 
the analogous equation for Uj. Therefore, we want to apply the classical result by Kruzkov [201 
Theorem 1], see also EH Theorem 2.1], to each equation in dSU, setting iteratively for i = 1 ,..., n 

f{t,x,u) = ipi{t,x,u,{'d*w){t,x)') and F{t,x,u) = (t,x,u,{i!) * w){t,x)) . 

To this aim, we check that the assumption (HI*) in EH Theorem 2.1], see also EOl Theorem 1], 
is satisfied. 

(HI*) / G C°(Dt;R^) holds by (if) and (-d), since w G C° (/; Li(R^;R")). 

F G C°(Dt;R) holds by ($) and {-&), since w G C° (/; (R^; R")). 

/ has continuous derivatives duf, du grad^, /, grad^ /, by (ip) and (i9). 

F has continuous derivatives and grad„, F by (€>) and (i?). 
d,/GL-(D^;R) by (cp). 

(R-div,/) G L-(0^;R) by (^p) and ($). 
i9„(F - diVa;/) G L°°(D^;R) by (ip) and ($). 

Therefore, problem m admits a solution u G L°° (/; LJqj,(R^; R”)). 
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(5.4) 


3: Total Variation Estimate. We want to apply [TUI Theorem 2.5] as refined in (UTJ Theo¬ 
rem 2.2]. To this aim, we verify (H2*) in [2T1 § 2]. 

(H2*) grad,,5„/ G by (y.) and (i9). 

duF G since 9„F = dui>fi and since (4>) holds. 

Ii /k" ~/)(^i OIIlooj-^j^.rn) dx dt < -boo: indeed, note that the inequality 

||grad,^F(t,a:, •)||Loo(,^j^.RAr) < X{t,x,U) holds by ($). Moreover, 

divx fit,x,u) = diyxy:>i{t,x,u,{'d*w){t,x)) 

+ grad^ ifi it, x,u, {d * w){t, x)) div^,(d * w)it, x) 

and passing to the gradient 

gradj, div:j,/(t,a:,w) 

= gradj, diva, tpi it,x,u,id * w)it,x)) 

+ gradyi diva; ipi (t,x,u,i'd * w)it,x)) grad^,(d * w)it, x) 

-b grad^, grad^ it,x,u,i'd * w)it,x)') div^, fd * w)it,x) 

+ grad^ifi it, x,u, it} * w)it,x)) grad,j,('d * rc)(t, x) diva;(^? + w)(t, x) 

-bgrad^ (fii [t,x,u, (i? + w)it,x)) grad,^, diva;(d + w)it,x), 
so that, using the standard properties of the convolution and ((p) 

llgradj, diva:/(t, X, u)|| < (1-b 3 Cc/i?T-b C;/T^) A(t, x, 17) 

< Cu il +RT + R^T^) Xit,x,U) 

and hence, using (<p), ($) and (15.3L 

I I Irn f)it,x, •)||l~(Wj,;R«) da;dt 

< (l|g^ad,,F(t,x,-)||Loo(i^„.RN) + ||gi'ad,,div,,/(t,x,-)||Loo(j^„.RN)) dxdt 

~ 11 it + RT + R^T^)\it,x,U)^dxdt 

= Cu il+RT + R^T^) AiT,U). 

To apply m Theorem 2.5], with reference to [UTJ (2.6)] compute first 
gvad^ duf it,x,u) = gTad,,du(piit,x,u,i'd*w)it,x)) 

+ grad^ duifi {t,x,u,i'd *w)it,x)) grad,^,(d *w)it, x) 

so that 

llgradj. 9„/llLc»(7xR'^xWt/;R'^x^) — + Cu CRT < Cu (1 + CRT) 

and 

«o = i‘^N + l)\\grad^duf\\^ “(7xRJVxW[/;RJVxN) 

xUu-,R) 

< i2N + l)Cuil + CRT) + Cu 

< CCuil + RT). 

Denoting Wn = /J^^^(cos d)^ dz?, use (15.5|) . (15.61) to obtain, for all t G I, 


(5.5) 


(5.6) 


TV(u,(t)) < TYiu)e<*+NWN [e^oit-r)f ||grad,(F - div,/)( t, x, •)! 

JI JR^ 

< (TYiu) + NWN [ [ ||grad,^(E-diV:j,/)(r,x,-)| 

\ Jl JR" 


L“>(Wi7;RN) 


dx dr 




dx dr 


1 e'"S* 


< (iG-bC'C[/(l-bi?T-bi?^r2) A(r,C/)) e^Cu(i+RT)T _ 


(5.7) 
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4: Continuity in Time We use [lU Corollary 2.4]. To this aim, verify first that f,F 

satisfy (H3*). 


(H3*) duf e already verified in (HI*). 

duF e L“(ng;R), already verified in (H2*). 

Ii /r" IK-^ “ ■)|Il“(W£/-R)^ 1I5-4I1 - (tp) and ($) to obtain the 

bound 


I Jr^ 


< 



Jl JRJV 

< A(r, u) + 




|diVx/(t,X,-)||L=o(;^^.I 



diVa; if (t, X, w)(t, x)) 

Jl JR!^ 

+ grad,i(p{t,x,-,i^*w)it,x)) 

Ji Jr!^ 

X IldiVa; ^?|Il~(7xR";R") II^IIc“(/;L1(k'' ; 

< A(r, u) + A(r, u) (i + Ijdiv, 

= C{IfRT)K{T,U). 


l°°(Wi7;R 
dx dt 


dx dt 


ll'“^llcO(/;Li(R";R))^ 


(5.8) 


Repeating the same computations on the time interval between s and t, by (y>), dH (2.8)], (15.7p 
and (|5.8I1 . for all t,s £ I, 


\ui{t) - 71^(5)11,1 


L1(RJV;R) 


< 


[ f ||(F-diva,/)(r,x,-)| 

Js Jr’^ 


, dx 


'L°°(Uu-, 

< C{l + RT)\K{t,U)-K{s,U)\ 

+C\t - s| [a: + C{1 FRT + R^T^)A(T, U) 


(5.9) 


+ \t-s\ xR^xWcJ.R) SUpTV(ui(T)) 

tGI 


^CCuil+RT)T 


(5.10) 


proving the uniform L^-continuity in time of the map t —>■ Ui{t), where u = F(w, u). 


5: L°° Bound Passing to the limit e —>■ 0 in the classical estimate [101 Formula (4.6)], we have 
that, using [101 Formulae (4.1), (4.2) and 4) in § 4], |■u(^,x)| < {Mo + CoT)e‘^^'^, where 


II^(^)|Il”(R'^;R) 

< R- 

diva,/(-,-,0) - F(-,-,0 )||lo„(^^jj„. 

diva, ipi {■,-,0,{'d*w){-,-)) 


Co = 


< 


L“(7xR'V;R) 

grad^ ifi (•, •, 0, (t? * w){-,-)) diva:(d + w){-,-) 


L“(7xRJV;i 


L“(7xR";: 


Cl = 


< 


$(-,-,0,(i?*'u;)(-,-)) 

< Cu + Cu CRT + Cjj 

< CCuil+RT). 

sup (-5„diva:/(t,x,?7) + 5uF(t,x,'u)) 


IxR^xUu 
dudivx ‘Pi (•,-,-,(i?*7y)(-,-)) 


L°°(7xR« xUu;'. 


14 






























Therefore, 


du grad^ <fi (r? * w)(-, •)) diva;(r? * w){-, •) 


L°°(IxR'^ xUu-,i 


L“(/xR'^xW[7;IB 


< Cu + Cu CRT + Cif 

< CCu{l + RT). 


\\uh^^,^^..^j,^<{R + CCu{l + RT)T)exp{CCu{l + RT)T) . (5.11) 


6: T is Well Defined Apply (I5.9I1 ~ (I5.10I) with s = 0, obtaining that for alH G / 


< 


(7(1 + RT)A{T) + CT K + C{1 + RT + R^T‘^)A{T) 




This inequality, together with (I5.11L ensures that if T is sufficiently small, u{t) = (T{w,u)) (t) G 
Bi^i (u, R) for all t G I. This estimate, together with what was proved at 2 and 4, ensures that 
T{w,u) G C° (/; L^(R^;R"-)) for any u G 


7: T is a Contraction. Here we use the stability result [H Theorem 2.6] as refined in [21] 
Theorem 2.5]. To this aim, for w',w" G C° (/; L^(R'^; R")), call f, f”, F', F" the corresponding 
fluxes and sources. We first verify that /' — /" and F' — F" satisfy (H3*). 

(H3*) duW - f") G L°°(flg;R^) is proved as in (HI*). 
du{F' - F") G L“(fl^;K) is proved as in (H2*). 

Using ($) and (|5.4p . 


/ 


|(A'-F")-div,(/'-/")(t,x,-)| 




dx dt 


(5.12) 


< 


< 


I JR«‘ 


\{F'-F"){t,x,-)\ 


L“(Wt7;: 


dx dt 


/ Jr" 


|div,(/' - •)||loo(j.^„.r) dxdt 


[ [ {t,x,-,{'d * w'){t,x)) - {t,x,-,{'d * w"){t,x)) 

JI Jr" 


L“(Wi7; 


dx dt 


div. 


I Jr" 


< 11$ 


wi.- 

Il</’llw2,c 


I Jr" 


ipi {t, x, •, (d * w'){t, x)) - (fit (t, X, ■,('&* w"){t, x)) 
grad^ (fi {t, X, ■,{'&* w')(t, x)) diva;(z9 * w'){t, x) 

- grad^ ipi (t, X, •, (r? * x)) div:j,(i? * w"){t, x) 

= (/xR"xWt/xW-;R")y^^^ -w")){t,x) dx dt 

■(/xR"xW( 7 XW™;R'“x;v) / / -w")){t,x) dxdt 

JI Jr" 

grad^ ipi (t, X, •, {-d * w'){t, x)) 


dx dt 


L“(Wc7;i 


dx dt 


L“(Wt/; 


- grad^ ipi (t, X, •, (-d * w”){t, x)) 


diva;('d * w'){t, x) 


dx dt 
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I Jr" 


grad^ V?i {t,x,-,{d*w”){t,x)) 

[div 2 ;(i? * w'){t,x) — div 2 :(i? * w"){t,x)\ 


da; dt 


L°°(Wt/;R) 

— II‘^’IIwi.“(7xR"xW[/xW5?;R")II^IIl“(R";R™x")||'^ II C0(7;Li(R";R"))'^ 

+ II‘/’IIw2.“(7xR"xW( 7XWJ}‘;R”X") ll^ll L“(R";R™xn) ||ui W || cO(7;L’-(R";R"))'^ 

+ ll‘/’llw2.~(7xR"xW[7XWJf;R"xiV) 

+ ll</’llwi.~(7xR"xW(7xW5‘;R"xiV)^^^ {dw^ i} * {w' - w")) {t, x) da; dt 

< CC'ar(l + i?)||77;'-77;"||^0(,^L,(^„^„„)). (5.13) 

Recall the following quantities from [211 (2-10)] and use (15.411 : 

K = ll^u-P" ||lcx=(7xrNxR;R) + ll^Ji ~ / ) ||l°°(7xR"xR;R) (5-14) 

— I|4’IIwi.“(R+xR"x77(/xWJ1*;R"xN) 

9„div:r (v’i * w')i-,-)) - <Pi {■,■,■,{■& * w")(-,-))) 

^ ' L“(7xR"x77t7;R) 

du grad^ (/Ji (•,•,•, (t? * w')(t)) diva;(7? *w'){-,-) 

-du grad^ ipi (•, •, •, {d * w"){-, ■)) div^(7l * w"){-, ■) 

L°°(7xR"x«t/;R) 

^ l|4’|lwi.“(R+xR"x77c/x77^;R’*xiV) + 2||;/?||-wr2,oo (r+^R" xWy xW^;R"xN) 

+2||V^|Iw2.“(R+xR" x77tf x775‘;R”XiV)l|diVa; ^Uloo (RN.Rn) ||?7; || cO(/;Li (R";R")) 

+ ll‘t^llw2.<x>(R+xR"x77(/x77^;R"xAf)||diVa; ^Uloo (RAr.R„) ||77; — W || co(/;Li (R";R")) 

< ^Cu + dCuC RT 

< CCu{l + RT). (5.15) 

M = ||9„/"|| 

L“(7xR" xR;R) 

duVi (•,-,-,(d*77;")(-,-)) 


L“(7xR" xWf/;R) 


< 


Iw2.°<=(r+xR" xWt/xW™;R"xiV) 


< Cu ■ 

By [5TJ Remark 2.8], (15.61) and (15.151) 

/SI * i /^ * i 

^ ° ~ ^ i ^ gC'Ct/(l+7iT)t 

^ ^ 

SO that we can prepare the bound 

\\duif —f )||L~(7xR"x77g) 

= du(fi (•, •, •, (d * «;')(•, •)) - duV>i (•,•,•, (d * M^")(') ■)) 

^ ll¥’*llw2’“(7xR"x77t/x77g?;R") ll^llL”(R";R"xm) lln; —W 

— C'C'gr IIit; —w ||co(7;L1(r";R)) ’ 


(5.16) 


L“(7xR" xUu) 

f I I 

llc“(7;Li(R";R)) 


(5.17) 
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and we can finally pass to the key estimate provided by pn Theorem 2.5] using (I5.16I1 . (I5.17I1 . 
(1^ . (1051) and (15:^ 


< 




K '*' _ lx* ’I' 

n ^ 


(5.18) 


+NWn I -;-:-||grad^(i^'-diV;r/')(r,a:,-)||Loo(2.,„.RiV)dxdT 


Kn — K 


x||5m(/ / )|lL“(7xR«‘xWt/) 


f [ {{F' - F") - div,(/' - /")) (t, X, •) 

Jo JR"‘ 


+ / e 
/o 

< KCCu 


L“(Wi7;R") 


dx dr 


C0(/;L1(R«‘;R)) 

TC + i?r + a(T, U) C Cc/ II w' - w 

/'II 


IC°(7;L1( 


< 


+ ^CC.(1+77T)^^^ T’ (1 + ^) Ih' - IIco(7;LM«-«)) 

C'C[/T(i + /s:+ (1 + i?T + i?2T2) A(r, U) + R^ ^cCuii+RT) 11^' 


— W 


ICO(/;Li( 




which shows that there exists a positive T*, such that the map 


Tu 


C°{[0,n];B^.{u,R,U)) 
w 


^ C°([0,T,];Bli(u,E,C/)) 

—^ T{w,u) 


is a contraction, for any u G i?LinBv(w, i?, U,K). 


8: The Fixed Point of T Is the Unique Solntion to (II.IL The fixed point of 7a solves (11.11) 
by Definition 12.11 and from (15.21) . On the other hand, any solution to (11.11) in the sense of Defini¬ 
tion O is a fixed point of 7a, proving also uniqueness. 


9: Continuous Dependence on the Initial Datum. Note first that T is L^-Lipschitz con¬ 
tinuous in its second argument. Indeed, applying again m Theorem 2.5], we have: 


\T{'W,u) T(w,U )||co(7;L1(R";R'‘)) — ® ||m U ||l1(r« 

< - u" 


lLl(R'V;i 


(5.20) 


By [3 Theorem 2.7], the L^-Lipschitz continuous dependence on the fixed point of 7a from u 
follows. □ 


Proof of Proposition 13.11 To prove that (13.81) fits into the class (HD), simply observe that 

* u = grad,j, (ry * (?7i + 772)). (5.21) 

The regularity required in (ip) and ($) is immediate, the cutoff function Fg being useful in 
bounding the terms grad^ p and grad^ p. The various estimates follow from (13.91) and from 
the fact that the map {xi,X 2 ) —>■ {xi,X 2 )/+ x 2 '^ is bounded, with all first and second 
derivatives also bounded. □ 


Proof of Proposition l47ll Observe first that (14.71) . the assumption e > e and (14.21) ensure that 
the flow in the convective part of (l+HI) points inward all along the boundary of B. Therefore, if 
there is a solution to (l+Hl) . its support is contained in B for all times. To apply Theorem 12.21 we 
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introduce a function s G C^(R^; R) such that s{x) = 1 for all x £ B. Then, note that (14.61) belongs 
to the class dLH. Indeed, similarly to (I5.21I1 . set 


N = 2 
^=1 

TO = 2 dx^ri{x) dx^r]{x) 

u = p 


^{t,x,u,A)=u(v{x)-^-^i^^^^^y{x) ^5 22 ) 

$(t, X, u, A) = x) - 4'o„t(<, x, u ), 


The invariance of B proved above ensures that the function s has no effect whatsoever on the 
dynamics described by (IMl)- Therefore, with the given initial datum (as well as with any other 
initial datum supported in ,B), any solution to (I1.1I1 ~ (I5.22I1 also solves (14.6L and viceversa. The 
estimates required in (ip) and (#) now immediately follow. □ 
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